t 


REPORT  DOCUMENTATION  PAGE 

form  Adprovad 

OM§  No.  0704^73M 

of  m#ontiftn<yi.  09%  mtMm  mom^m  or  wm  «odct  of 

- - TTUin  nr  -mn--t  to  '*«Mdu«rton  Vtfvicot.  Oirocwto  for  mformooow  QoTooom  o<%<  dra.  lill  ioftofw 

tnoOftKt  Of  lod  iMOoot.  Amctiod^rotdct (07044 OC  WWX 

t.  ACENCf  USE  QNtt  (tMvm  Oittuc) 

2.  RERORT  DATE 

01/04/02 

V/VuWmL 

4.  HTU  ANO  SUITTOE 

Electromagnetic  Scattering 

S.  7UN04NG  NUMRERS 

' 

F49620-99-1-0010 

E.  AUTMOmS) 

Oscar  P.  Bruno 

7,  PCAAOEMiNO  OAGANUATION  NAME<S)  ANO  AOORESS(ES) 

Applied  &  Computational  Mathematics 

California  Institute  of  Technology 

Pasadena,  CA  91125 

S.  PERECRMING  ORGANIZATION 
REfCRT  NUMftER 

SPONSOWNa/MONITOWNG  AGENCY  NAM£(S)  ANO  AOORESS<ES) 

USAF,  AFRL 

AF  OFFICE  OF  SCIENTIFIC  RESEARCH 
,  801  N.  RANDOLPH  ST.  ROOM  732 

ARLINGTON,  VA  22203-1977 

10.  SPONSORING/  MONITORING 

AGENCT  REPORT  NUMER 

11.  SUmXMINTAllY  NOTES 


12a.  tMSTWWnOM/AVAItAWUTY  STATEMENT 

Unlimited.  Approved  for  public  release. 


1 12ft.  OtSTWIUTION  coot 


IX  AtSTMCT  (MsjamumJOOwonixJ  ■,  t  ■  ‘  ' 

*  Application  of  the  method  of  boundary  variations  to  the  calculation  or  elect  ; 

magnetic  scattering  from  the  ocean.  _  _  ' 

*  Introduction  of  new,  high  frequency  approximations  (series  expansions  in  powers  : 
of  the  wavelength)  for  use  in  conjunction  with  boundary  variations  in  ocean 

applications.  ^  .  i 

*  Application  of  the  method  of  boundary  variations  in  the  solution  of  eigenvalue 

probiems.  _  '  ;  \ 

*  Development  of  fast  high-order  integral  methods  for  solution  of  problems  of 
electromagnetic  scattering  involving  surfaces  and  volumes  in  three  dimensional  o;) 
space. 


20020130  309 


14.  SURJECT  TERMS 

IS.  NUMSER  OP  PAGES  ^ 

1E.PRIUCOOE 

17.  SECURITY  OASSIPKATION 

OP  REPORT 

11.  SECURITY  CLASSIPICATION 

OP  THIS  PAGE 

1 S.  SIGJWr  JlSSiWCATlOfl 

:0^  Aisiiucr 

20.  UNRTADON  OP  AESTRACT 

Unlimited 

NSN  7540^1-280-5500 


Standard  form  298  (Rav.  2*89) 


GENERAL  INSTBUCnONS  FOB  COMPLETING  SP  298 


Th*  R«oori  Documentation  Page  (ROP)  iJ  used  in  announcing  and  cataloging  reports,  it  is  important 
that  this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  title  page. 
Instructions  for  filling  in  each  blodt  of  the  form  follow,  it  is  important  to  stay  wit/un  the  //nej  to  meet 

optfca/ scanning  reguireflieno. 


Block  1.  Aoencv  Use  Only  (Leave  b/ankj. 

Block!-  Report  Date.  Full  publication  date 
including  day,  month,  and  year,  if  available  (e.g.  1 
Jan  88).  Must  cite  at  least  the  year. 


Block  12a.  Oistnbution/Availabilitv  Statement. 
Denotes  public  availability  or  limitations.  Gte  any 
availability  to  the  public.  Enter  additional 
limitations  or  special  markings  in  ail  capitals  (e.g. 
NOFORN,  REL,  ITAR). 


Blocks.  Type  of  Report  and  Oates  Covered. 
Statewhether  report  is  interim,  final,  etc.  if 
applicable,  enter  inclusive  report  dates  (e.g.  10 
Jun87-30Jun88). 


Block  ♦.  Title  and  Subtitle.  A  title  is  taken  from 
the  part  of  the  report  that  provides  the  most 
meaningful  and  complete  information.  When  a 
report  is  prepared  in  more  than  one  volume, 
repeat  the  primary  title,  add  volume  number,  and 

indude  subtitle  for  the  specific  volume.  On 
dassified  documents  enter  the  title  d  issification 
in  parentheses. 


Blocks,  Funding  Numbers.  To  induce  contract 
i.nd  grant  numben;  may  indude  program 
element  number<s),  project  number($),  task 
numberfs),  and  work  unit  numbef(s).  Use  the 
following  labels: 


C  -  Contract 
G  •  Grant 
R€  -  Program 
Element 


PH  -  Project 
TA  -  Task 
WU  -  Work  Unit 

Accession  No. 


Blocks.  Authorfs).  Name<s)  of  person(s) 
responsible  for  writing  the  report,  performing 
the  research,  or  credited  with  the  content  of  the 
report.  If  editor  or  compiler,  this  should  follow 
the  name($). 


Block?.  Performing  Organization  Name(s)  and 
Address! es).  Self-explanatory. 

Block  1.  Performing  Oroanitation  Report 
Number.  Enter  the  unique  alphanumeric  report 
numberfs)  assigned  by  the  organization 
performing  the  report. 

Block?.  Soonsorino/Moni  tori  no  Agency  Name(s) 
and  Addressles).  Self-explanatory. 


Block  10.  Soonsorino/Monitoring  Agency 
Report  Number.  (If  known) 

Block  11.  Supplementary  Notes.  Enter 
information  not  included  elsewhere  such  as: 
Prepared  in  cooperation  with...; '  rans.  of...;  To  be 
published  in....  When  a  report  is  revised,  include 
a  statement  whether  thu  new  report  supersedes 
or  supplements  the  older  report. 


DOO  -  See  DoOO  5230.24, 'Distribution 
Statements  on  Technical 
Documents.* 

DOE  •  See  authorities. 

NASA  -  See  Handbook  NH8  2200.2. 

NTIS  -  Leave  blank. 


Block  12b.  Distribution  Code. 

DOO  •  Leave  blank. 

DOE  •  Enter  DOE  distribution  categories 
from  the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports. 

NASA  •  Leave  blank. 

NTIS  -  Leave  blank. 

Block  13.  Abstract.  Indude  a  brief  fJUax/mum 
200  words)  factual  summary  of  the  most 
significant  information  contained  in  the  report. 


Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 


Block  15.  Number  of  Pages.  Enter  the  total 
number  of  pages. 

Block  16,  Price  Code.  Enter  appropriate  price 
code  (NTIS  on/y). 

Blocks  17.- 19.  Security  Oassifications.  Self- 
explanatory.  Enter  U.S.  Security  Qassification  in 
accordance  with  U.S.  Security  Regulations  (i.e., 
UNCLASSIFIED).  If  form  contains  classified 
information,  stamp  dassification  on  the  top  and 
bottom  of  the  page. 


Block  20.  Limitation  of  Abstract.  This  block  must 
be  completed  to  assign  a  limitation  to  the 
abstract.  Enter  either  UL  (unlimited)  or  SAR  (same 
as  report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited.  If  Wank,  the  abstract 
is  assumed  to  be  unlimited. 


Staneard  Form  iSe  la<x  (R«v.  2  «9) 


Submitted  via  email  Jan  4,  2002 

Dear  Arje, 

Below  please  find  the  final  report  for  my  AFOSR  award 
#F49620-99-l-0010 . 

Would  you  please  acknowledge  your  receipt? 

Best  regards, 

Oscar  Bruno 

Electromagnetic  Scattering 
Grant  #F49620-99-l-0010 
Final  Report 
Oscar  P.  Bruno 
January  2002 


OBJECTIVES 

*  Application  of  the  method  of  boundary  variations  to 
the  calculation  of  electromagnetic  scattering  from  the 
ocean 

*  Introduction  of  new,  high  frequency  approximations 
(series  expansions  in  powers  of  the  wavelength)  for 
use  in  conjunction  with  boundary  variations  in  oceanic 
applications 

*  Application  of  the  method  of  boundary  variations  in 
the  solution  of  eigenvalue  problems 

*  Development  of  fast  high-order  integral  methods  for 
solution  of  problems  of  electromagnetic  scattering 
involving  surfaces  and  volumes  in  three  dimensional 
space 

RECENT  ADDITIONS  TO  THE  RESEARCH  PROGRAM 

*  Development  of  a  perturbative  technique  for  the  normal  mode  problems 
associated  with  lasing  cavities  (extension  of  work  on  eigenvalue 
problems) . 

*  Development  of  a  perturbative  technique  for  the  numerical  solution 
of  the  non-linear  Schoedinger  equation  for  the  propagation  of  pulses 
in  optical  fibers. 

*  Solution  of  inverse  problems  related  with  Optical  Coherence  Tomography 
(Initial  stages,  in  collaboration  with  graduate  student  Julian  Chaubell. 
This  work  is  being  supported  in  part  by  a  grant  from  the  TRW  foundation.) 
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REPORT 

BOUNDARY  VARIATIONS /HIGH  FREQUENCY  APPROXIMATION  FOR  ROUGH  SURFACE 
SCATTERING:  MULTI-SCALE  SOLVER 

Our  approach  to  the  solution  of  rough-surface  problems  containing 
multiple  scales  (including  oceanic  scattering)  is  based  on  a 
combination  of  the  ""Method  of  variation  of  boundaries’’,  developed 
earlier  in  collaboration  with  F.  Reitich,  and  a  high-order 
high-frequency  approximation,  developed  recently  in  collaboration  with 
A.  Sei.  Our  last  year’s  efforts  allowed  us  to  fully  develop  and 
demonstrate  the  capabilities  of  the  combined  high-order  high-frequency 
solver  with  the  boundary  perturbation  approach.  The  complete 
multi-scale  algorithm,  including  boundary  variations  and 
high-frequency  solutions,  has  been  implemented  in  a  Fortran 
code.  Applications  to  various  configurations  of  interest  are  currently 
underway. 

Additional  work  in  the  area  of  perturbative  methods  concerns 
normal-modes  in  cavities.  This  work,  in  collaboration  with 
undergraduate  student  Arnaud  Durand,  generalizes  the  perturbative 
method  for  eigenvalue  problems  of  Bruno  and  Reitich  to  the  physically 
relevant  normal  modes,  that  is,  eigenvalue  problems  in  all  of  space, 
(Mr.  Durand  spent  the  summer  of  2001  visiting  our  group  from 
Cachan,  Paris.) 


FAST  INTEGRAL  EQUATION  SOLVERS:  VOLUMETRIC  SCATTERERS 


The  postdoctoral  student  McKay  Hyde  is  working  under  my  direction  in 
the  development  of  fast,  high  order  solvers  for  problems  of  volumetric 
scattering.  We  currently  have  a  fully  tested  version  of  our  2d  code, 
as  well  as  a  three  dimensional  volumetric  solver  for  smooth 
scatterers.  This  three  dimensional  solver  has  been  implemented  so  as 
to  run  in  a  massively  parallel  environment.  Current  work  in 
collaboration  with  F.  Reitich,  which  is  nearing  completion,  will 
provide  the  elements  necessary  to  extend  this  solver  to  discontinuous 
scatterers.  These  solvers  will  be  tested  against  experimental  results 
given  in  the  thesis  of  Jeffery  M.  Parks,  Captain  USAF  -  kindly 
provided  to  us  by  Andy  Greenwood  from  AFRL. 

FAST  INTEGRAL  EQUATION  SOLVERS:  SURFACE  SCATTERERS 

Over  the  last  year  L.  Kunyansky  and  I  a  class  of  *high-order  methods 
for  singular  surfaces*,  that  is,  an  algorithm  which  exhibits  high 
order  convergence  even  for  surfaces  containing  corners,  edges,  etc.  We 
have  demonstrated  high-order  convergence  in  the  well  known  ogive 


geometry  of  EMCC;  indeed,  we  have  computed  scattering  from  a  20*lambda 
ogive  with  several  digits  of  accuracy  in  a  three-hour  computation  on  a 
400  Mhz  PC  (with  1Gb  of  memory) .  To  the  best  of  our  knowledge  such 
accuracies  are  absolutely  unprecedented  for  a  singular  geometry  of 
this  type,  of  a  size  as  large  as  the  one  we  considered;  it  is 
specially  encouraging  to  find  the  required  computation  time  was  very 
small.  A  somewhat  particular  code  was  produced  for  singular  surfaces, 
which  showed  that  our  strategies  for  edge  problems  are  sound.  Our 
p2:esent  work  seeks  to  produce  a  general  code  which  can  automatically 
incorporate  CAD  data,  and  can  resolve  geometric  singularities  without 
human  intervention.  This  research  is  being  conducted  in  contact  with 
Lockheed  Martin  and  the  Jet  Propulsion  Laboratory.  These  solvers  will 
be  tested  against  experimental  results  which  will  be  provided  to  us  by 
Andy  Greenwood  from  AFRL,  as  well  as  experimental  data  on  models  for 
the  Mars  lander  provided  by  Vaughn  Cable,  from  JPL.  In  a  final 
related  effort,  we  are  developing,  in  collaboration  with  undergraduate 
student  Stephane  Lintner,  integral  equations  suitable  for  the 
treatment  of  thin  plates  (such  as  a  disc) .  These  integral  equations 
should  play  an  important  role  in  the  treatment  of 

air-vehicles.  (Mr.  Lintner  spent  the  summer  of  2001  visiting  our  group 
from  Cachan,  Paris,  and  he  plans  to  join  our  group  as  a  graduate 
student  in  2002.) 

OPTICAL  COHERENCE  TOMOGRAPHY 

Optical  coherence  tomography  is  a  non-invasive  imaging  technique  based 
on  the  use  of  light  sources  exhibiting  a  low  degree  of  coherence.  Low 
coherence  interferometric  microscopes  have  been  successful  in 
producing  internal  images  of  thin  pieces  of  biological  tissue; 
typically  samples  of  the  order  of  1mm  in  depth  have  been  imaged,  with 
a  resolution  of  the  order  of  10  to  20  microns  in  some  portions  of  the 
sample.  Such  images  have  been  produced  through  renderings  of  the 
intensities  of  certain  interference  fringes  as  functions  of  the 
position  of  the  light -focus  within  the  sample;  quite  generally, 
limited  post-processing  of  this  data  has  been  used.  In  our  work,  in 
collaboration  with  graduate  student  J.  Chaubell  (under  partial  support 
of  TRW)  we  address,  in  a  mathematically  rigorous  manner,  the  inverse 
(Maxwell)  problem  of  producing  the  actual  values  of  the  refractive 
index  within  the  sample  from  given  low-coherence  interferometric  data. 
Once  obtained,  such  a  map  of  the  refractive  index  variations  is  useful 
in  a  variety  of  ways;  in  particular,  a  straightforward  display  of  this 
map  yields  an  image  of  the  internal  structure  of  the  sample.  Our  work 
over  the  last  year  has  produced,  for  the  first  time,  a  well  posed 
inverse  solver  based  on  optical  coherence .  We  expect  important 
developments  may  follow  from  these  contributions. 

PERTURBATIVE  SOLUTION  OF  THE  NON-LINEAR  SCHOEDINGER  EQUATION  FOR 
MODELLING  OF  OPTICAL  FIBERS 

In  collaboration  with  former  postdoctoral  associate  D.  Amundsen  (now 
assistant  professor  at  Carleton's  mathematics  department  (Canada))  we 
developed  a  perturbative  technique  for  the  solution  of  the  fully 
nonlinear  Schroedinger  equation.  This  work,  which  was  motivated  by 
contacts  with  Corning  (NY)  and  Fiberspace  (CA)  has  produced  a  highly 
efficient  solver  -  which  can  speed  up  computations  by  as  much  as  one 
hundred  times.  This  is  a  significant  accomplishment,  since  some 
industrial  applications  require  repetitive  solution  of  these 


equations,  each  one  of  which  requires  of  the  order  of  10  hours.  This 
technique  should  be  useful  in  other  nonliner  evolution  problems,  and, 
we  believe,  it  generally  has  a  very  significant  potential  indeed.  The 
basis  of  the  method  is  to  utilize  analytic  continuation  *instead*  of 
finite  differences  in  the  numerical  method.  Thus,  much  longer  time 
steps  can  be  used  in  an  algorithm  which  is  essentially  explicit,  and 
at  the  same  time  unconditionally  stable. 

ADDITIONAL  EFFORTS 

In  addition  to  these  efforts  related  to  electromagnetic  scattering,  in 
the  last  year  this  PI  was  involved  in  a  variety  of  research  efforts  in 
the  general  area  of  materials  science,  including  work  on 
graphite- to-diamond  phase  transitions  induced  by  shock  waves  (in 
collaboration  with  postdoctoral  associate  Dimitri  Vaynblat,  and,  most 
recently,  with  Prof.  Ruben  Rosales  (MIT)  and  his  student  John 
Wheatherwax) ,  work  on  magnetoelastomers  (in  collaboration  with 
Prof.  Liliana  Borcea,  Appl.  Math.,  Rice  University);  and  solution  of 
free  boundary  problems  arising  from  finance  (in  collaboration  with 
graduate  student  Andy  Greenberg) . 
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